Load packages
suppressWarnings(suppressMessages({
library(Seurat)
library(magrittr)
library(imager)
library(EBImage)
library(STutility)
library(magrittr)
library(dplyr)
library(DT)
library(kableExtra)
}))
SVZ.Visium <- readRDS("SVZ.Visium/SVZ.Visium.rds")
Normalization
suppressWarnings(suppressMessages({SVZ.Visium <- SVZ.Visium %>%
SCTransform(verbose = F) %>%
RunPCA(verbose = F) %>%
RunUMAP(reduction = "pca", dims = 1:20,verbose=F)}))
Idents(SVZ.Visium)="labels"
SVZ.Visium <- RenameIdents(SVZ.Visium,"svze18"="E17.5","svzp2"="P2","svzp12"="P12","svzp22"="P22")
SVZ.Visium$labels <- Idents(SVZ.Visium)
SVZ.Visium@meta.data[["labels"]]=as.character(SVZ.Visium@meta.data[["labels"]])
Visualize only the spots
ST.FeaturePlot(SVZ.Visium, features = c("Sox2", "Pax6"), ncol = 2, grid.ncol = 2, pt.size = 4, show.sb = FALSE)
## Le chargement a nécessité le package : viridis

Visualize normalized gene expression in spots within sections
FeatureOverlay(SVZ.Visium, features = "Sox2",
sampleids = 1:4,
cols = c("light grey", "mistyrose", "darkred"),
pt.size = 3,
add.alpha = TRUE,
ncol = 4, show.sb = FALSE,
pt.alpha = 0.1,label.by = "labels")

Visualize normalized gene expression in spots with section
FeatureOverlay(SVZ.Visium, features = "Lrp2",
sampleids = 1:4,
cols = c("light grey", "mistyrose", "darkred"),
pt.size = 3,
add.alpha = TRUE,
ncol = 4, show.sb = FALSE,
pt.alpha = 0.1,label.by = "labels")

Other visualization parameters
FeatureOverlay(SVZ.Visium, features = "Sox2",
sampleids = 1:4,
cols = c("lightgray", "mistyrose", "red", "darkred", "black"),
pt.size = 1.5,
add.alpha = TRUE,
ncol = 4,label.by = "labels")
